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ABSTRACT 

A CJalerkin method for systems of PDE’s in circular geometries is presented with motivating 
problems being drawn from structural, acoustic and structural acoustic applications. Depend- 
ing upon the application under consideration, piecewise splines or Legendre polynomials are 
used when approximating the system dynamics with modifications included to incorporate 
the analytic solution decay near the coordinate singularity. This provides an efficient method 
which retains its accuracy throughout the* circular domain without degradation at the singu- 
larity. Because the problems under consideration are linear or weakly nonlinear with constant 
or piecewise constant coefficients, transform methods for the problems are not investigated. 
While the specific method is developed for the 2-D wave equation on a circular domain and the 
equation of transverse motion for a thin circular plate, examples demonstrating the extension of 
the techniques to a fully coupled structural acoustic system are used to illustrate the flexibility 
of the method when approximating the dynamics of more complex systems. 
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1 Introduction 


A common step in many structural, acoustic and fluid applications involves t he solution of par- 
tial differential (‘((nations (1M)F s) modeling the physics of a system having a circular geometry, 
hxamples of applications in which this is important include the determination and control of 
circular plate* and cylindrical shell dynamics, tlu* modeling and attenuation of noise* within a 
cylindrical cavity, the simulation of fluid flow and boundary layer growth in a pipe*, as well 
as innumerable other applications involving circular domains. Due to tin* complexity of the 
problems, however, analytic solutions usually cannot be obtained and one must numerically 
approximate the dynamics under consideration. 

from separation of variables and truncation of the resulting infinite series, approximations 
in the form of modal expansions involving Bessel components can occasionally be used to 
approximate RDF dynamics. However, modal expansions can l>e determined for only a limited 
number of applications and hence one must often employ moregenera.1 techniques such as finite 
difference, finite element, spline, or spectral expansions in order to approximate the dynamics 
of the system under consideration. 

In the category of spectral methods for circular geometries, (lalerkin. collocation and tan 
methods have been studied with the choice of method depending upon the problem being con- 
sidered. lo date*, much of this research has centered around the* simulation of fluid flow and 
boundary layer growth and in these cases, emphasis has usually been placed on collocation dm* 
its success in handling complex boundary conditions, variable* coe*fficie*nts, and strong nonlin- 
e'arities [l(>. 1 7]. (lab* rkin methods for flows on spheres are* discussed in [lb], but this is done* 
primarily in the* context of Fourier expansions involving the* strong form of the* modeling flow 
equations with only a brief discussion concerning Legendre base's be*ing included. As noted in 
this latfeT reference, care* must lx* ex<*rcised whe*n applying any of the* previously mentione'd 
techniques ((lalerkin, collocation or tan) to problems with coordinate singularities since* the* 
incorrect application of pole* conditions can significantly degrade the* aecuracy of t.lu* me'thod 
as we'll as introduce* strong nunu'rieal instabilities. This refeivnce also includes a general com- 
parison bet ween the* results obta.in<*d with surface* harmonics (eigeuifunctions of tlx* Laplacian), 
modilied Robert functions and Fourier se*ries using collocation and (Iale*rkin methods in the* 
presence of the* coordinate* singularity. The* use* of the* modifie'd Robert futetions and tech- 
niques lor improving their conditioning and e*mph)ying fast transforms is furthe*r addressee! in 
[10]. finally, general ove*rvie*ws conce*rning the* application of spectral netluxls to prol)le*ms 
wit h coorelinate* singularit ies can be* founel in [II, 12]. We point out t hat the* ennphasis in tlu*se* 
re*f(*re*nce*s (enters primarily around the* de*V(*lopme*nt of nume*rieal methods for flows of various 
types and hence the* techniques are* often ge*ared toward t he approximation of st rongly nonline*ar 
operators which requires transleirm me'thexls to facilitate* e*flicie*nt implementation. While* these* 

1 eclmiques will clearly week for line*ar problenns involving circular domains, the *v e)fte*n are* not 
optimal nor the* be*st, choiee for the* applications unde*r eensidenat ion. 

In this paper, we* presemt a (!ale*rkin nethexl for liiear or vve*a.kly nonlinear probl(*ms having 
circular or cylindrical domains using pie*re*wise* spline* and spectral bases. Two aieas from which 
we* will draw <*xample*s are* structural elynainics and acoustics (a se*arch of the lit erat me reveals 
fe*vv numerical techniques for these* type’s of problems on circular domains). Linearization in 
the* first, ease* is olten justified when d<*aling with small amplitude* vibrations while* linearized 



acoustic equations are often employed when considering acoustic fields having sound pressure 
levels less than 150 dB (which is the case in a large number of acoustic and structural acoustic 
applications). While weak nonlinearities of the type considered in [7] can also be efficiently 
implemented, the discussion here will concentrate on the linear case. Moreover, the coefficients 
in the problems of interest will be taken to be constant or piecewise' constant in nature (this is 
certainly a reasonable assumption for many acoustic problems and a valid condition in a large 
class of structural applications). Hence we will not address the use of transform methods for 
evaluating variable coefficients and nonlinear components. 

Throughout the discussion, the modeling equations will be approximated in the weak or 
variational form in accordance with energy formulations of the problems. In the* structural 
application, this is done to reduce smoothness requirements on approximating elements and 
to accommodate structural and material discontinuities as well as unbounded (discontinuous) 
input operators. The consideration of the acoustic problem in weak form proves to be useful 
when considering coupled structural acoustic systems. Moreover, in weak form, some physical 
boundary conditions (e.g., the hard wall conditions which we consider here) are natural which 
implies that no basis alterations are necessary when implementing the method. We point out, 
however, that if the application warrants, approximation of the acoustic dynamics using the 
strong form of the system equations can just as efficiently be accomplished using the techniques 
of the paper. Moreover, these techniques can be easily modified to accommodate a tail method 
if dealing with more complex essential boundary conditions (boundary conditions which must 
be explicitly satisfied in a variational formulation of the problem). 

In the second section of this paper, a 2-1.) acoustic problem on a circular domain is con- 
sidered with a Fourier- Legendre basis being used in the ensuing approximations. In addition 
to its accuracy, the use of this basis proves to be advantageous in control applications since it, 
facilitates the maintenance of uniform stability margins under approximation [•!]. Moreover, 
the the use of the spectral basis facilitates approximation in a quotient space' which is often the 
natural state space in acoustic problems. 

The techniques are then applied to a structural application in Section 5 with the problem 
of approximating circular plate dynamics being used to motivate the analysis. In this case, 
cubic splines modified to satisfy boundary conditions are used in conjunction with a Fourier 
basis to obtain a suitable finite-dimensional approximation to the problem. The use of splines 
rather than a spectral basis is motivated by smoothness requirements as well as the ease with 
which they can be adapted to satisfy essential boundary conditions, and can easily be changed 
if warranted by the form of the modeling equations and boundary conditions. We emphasize 
that when developing a general technique that could be applied in both case's, rare was taken to 
treat the coordinate singularity in a manner that prevented degradation of accuracy using either 
basis (spectral or piecewise splines) as well as avoided the introduction of spurious numerical 
instabilities. 

The techniques from Sections 2 and 3 are then combined in Section 4 where the problem of 
approximating the dynamics of a fully coupled structural acoustic system is considered. 1 his 
demonstrates the flexibility of the method when studying a more complex coupled physical 
system and illustrates the manner through which the basic techniques can be extended to more 
general problems. The convergence properties of the method and its use in determining physical 
properties of the system are illustrated through a, set of examples. Finally, advantages of the 
method in control applications are discussed in the concluding remarks. 



2 The Wave Equation 


In this section we consider tin' 2-L) wave equation on a circular domain il of radius a with tin 1 
boundary denoted by F. For Neumann boundary constraints and initial conditions 0 O and <j>i, 
the equations of motion art* 

(fht = c 2 A<j) + <j(t,r,0) , (?\ 0) e n J > 0 , 

v^*r = 0 , {rj) 6 r .f >0 , 

( 2 . 1 ) 

0(0, r,0) = 0o(r, 0) , (?*, 0) G , 

0 t (O, r, 0) = <i>](i\0) , (;\0) G ^ 


where 0, in acoustics applications [14], d(*notes tilt* velocity potential and r is the speed ol 
sound. We note* that tin* Nruinann boundary conditions are chosen so as to be consistent with 
the hard wall conditions in the physical application described in Section 1 and can bo altered 
to fit the individual problems under investigation. 

To gain insight into t in* behavior of the solution and to motivate the form of the approx* 
imation, wo first solve the equations analytically through a standard separation of variables. 
For the homogeneous problem {(} — 0), we take <f)(t,r,6) — 7 ( / ) <t> ( 7% 0) to arrive at the two 
dimensional Ilelmholz equation 

A + 7 2 4> = 0 , (r, 0) G n , 

(2.2) 

v 4> . r — o , {i\ 0 )e r 


and tlu* rc*l at ion T" + u> 2 T = 0 , t > 0 . The separation constant here is —j 2 = — (~) 2 where u> 
is the circular frequency with units of radians/sec. To find 4y we separate variables once more. 
Letting 0) = ft(r)0(0), the expansion of Helmholz’s equation yields the expression 

1 d ( dH\ 1 <P& 2 n 

vRdr V dr ) r 2 0 dO 2 


which implies that R and 0 must satisfy the differential equations 

ft {a) = 0 , 


R = 0 




and 


0" + m 2 e = 0 

0(0) = 0(2tt) . 


(2.4) 


The general solution to (2.4) is periodic and is given by 0 ( 0 ) = Ait " n9 , m = 0, ±1, ±2, • • • . 
Moreover, the general solution of the Bessel equation (2.4) is R(r) = Ai-Lnil 1 ') + A, }V„ t (7 r ) 
where ./,,, and Y m are the m th Bessel functions of the first and second kind, respectively. As 
noted in [1], these functions display the asymptotic behavior 


2 m w! 


Y m (~) 


2 m ( m - 1)! 


(2.5) 
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for fixed m and c — > 0. Because Y m grows unboundedly at the origin, we take /t : $ = 0 in order 
to guarantee that the solution remains bounded at r = 0. The 4 de'cay properties ol J, n will lx 4 
incorporated when developing an ap])roxiination method for the problem. 

The eigenvalues 7 mn of (2.2) are tin'll determined by applying the boundary condition and 
solving for the zeros of the nonlinear equation 


dd ni ( *7 ^ ) 
dr 


= 0 . 


(2 .fi) 


The corresponding eigenvectors or modes are 0) = ,/ m (7 lor tn = 0. ± 1. ±2, * • •, 

n = 1,2, ‘b ■ * • which implies that standing waves have the form 


(f>mn{t,r\0) = [« n sin (7 mn ct) + h n cos( 7 m?i c/ .)] J m ('y mn r)e 


ituO 


A variation of parameters approach can then be used to extend this to tin* nonhomogeneous 
case. With suitable smoothness assumptions on tin* initial conditions and forcing term, it can 
be shown that the solution to (2.1 ) has the form 

0( LrJ))= £ jr <t>mn ( / )./,„( %nnV)r lm0 ( 2.7) 

m= — oo it — 1 


where the coefficients d> mn {t) depend on and the 4 time'-de'peMide'nt Fourier* coefficients of 

(j (we refer the* reader to any standard PL)E text for derails concerning the justification of tin' 
formal calculations outlined here). 

Having outlined the arguments leading to the analytic form of the solution, we now want to 
use those results to motivate' an appropriate CJalerkin technique for approximating the solution. 
To facilitate the use of this scheme in coupled structural acoustic systems (see, for example, 
the system discussed in Section 4), the weak form of the system equations will be considered. 
Moreover, to illustrate' the use of t he method when approximating the acoustic' potential (which 
is determined to within only a constant [14]), a basis suitable lor a quotient space will be chosen. 

To pose the system (2.1) in weak form, we take the 4 state* to lx 4 <f> in the space II = L 2 (il) 
where L 2 ($i) is the quotient space of L 2 over the constant functions. Flu 4 space' of test functions 
is taken to be the quotient space V — The' inner products for the two spaces are 4 taken 

to be 


(^a)//= / , {<!>,£) V = / V<£ • V£f/u? 

Ju c* . hi 

(the overbars here 4 denote' roni|)l('X conjugation as com])are'd with those 4 used above 4 to delineate 4 
the quotient space's). Energy considerations or integration by parts them yields the variational 

forni - . _ f — f - 

<l>tt£dtj + j V<7; ■ V£e/uj = j (j^dto (2.8) 


for all £ in V . 

To discretize, we 4 bc'gin with a Fourier e'xpansion in 0 which yields the 1 approximate' solution 


4> M (t,v,0)= T (j), n (l ,r)( im0 . 

m=-M 


(2.9) 


4 



We point out that the use of the complex Fourier expansion simplifies the following discussion 
both in describing the form of the approximate solution and the construction of the system 
matrices. However, when combining these wave results with those of the circular plate to 
yield an approximation scheme for the coupled system (see Section 4), it is easier to use a 
real Fourier expansion when performing the actual computations (this is due to the presence 
of the piezoceramic patches on the plate). The interchange between the two expansions is 
straightforward, and hence details concerning the implementation of the real Fourier scheme 
are left to the reader. 

Several possibilities exist for spline, Legendre or Chebyshev expansions of <f) m (t.r) both in 
a collocation and (Jalerkin setting. These include direct expansions which maintain the parity 
of the solution, refinements to incorporate the decay of the solution at the origin, and mapped 
expansions which use all the polynomials and yield better center resolution. These expansions 
must satisfy the condition 

<H M 

r\ 


at the origin which guarantees the uniqueness of the solution. This yields the requirement 

= 0 at r = 0 , 77i ^ 0 . (2.10) 


To guarantee differentiability at the origin, it is appropriate to requin' that the remaining 
component satisfies 


dr* 


= 0 


at r = 0 . 


( 2 . 11 ) 


As detailed in [6], one expansion of </> m (/,r) which satisfies these properties is 


where 


and 


<MM*) 


jY m 


?i=0 


771 = 


771 


|m| = 0, • • • , 5 

M = 6, * * • , M 


I7(r) 


f\(r) —1/3 , 777 — 0, n = 1 

P n ( r ) , otherwise . 


( 2 . 12 ) 


Here / J n (r) is the n ih Legendre polynomial which has been mapped to the interval (0,«). The 
term / ; i(/’) — 1/3 when 777, = 0,7? — 1 results from the orthogonality properties of the Legendre 
polynomials and arises when enforcing the condition f n r, 0, z)du: — 0 so as to guarantee 

that the functions art' suitable as a basis for the quotient space. The inclusion of the weight 
r l m l incorporates the decay of the analytic solution near the origin (see (2.5)) while ensuring 
its uniqueness at that point. Finally, we note that if N Legendre functions are used, the limit 
N m is givt'n by N m = N + 1 when |m| ^ 0 and N m = N when 777 = 0, which implies that 
j\4 = (2 M + 1 )(N m + 1) — 1 basis functions are used in the wave expansion. 

Summarizing, the approximate solution to (2.8) can be expressed as 


M N m 

7n = — 1 V/ n = 1 


EMOBf(r,0) 


k= 1 


(2.13) 



where P™ (r)e ,md . In comparing the form of this approximate solution with the 

analytic wave solution defined in (2.7), it can be seen that we have essentially replaced the Bessel 
components by weighted Legendre functions. This provides an approximate solution which is 
more flexible and is quicker to converge in a variety of applications than that obtained with a 
purely eigenfunction expansion (indeed, for more complex systems involving wave components, 
the eigenfunctions are unavailable and would have to be approximated before they could be 
used as a basis). 

To provide an approximating finite-dimensional space for the problem, we define the sub- 
space H M = span{B J k M } and take the product space for the first-order problem to be H M — 
H m x H M . The restriction of the infinite-dimensional system (2.8) to the space H M x H M 
then yields 

for £ in H M . The corresponding matrix system is 

' o i r d M {t) I r o 
0 J [ # M (t) J ” [ ~ K r 

' o j r ^(o) j _ r 

0 M# J [ ^(0) J _ [ g? 

where i J M (t ) = [(j>\(t), <f> 2 (t), • • • , <!>m{^)Y denotes the M X 1 vector containing the approximate 
state coefficients. The component matrices and vectors are given by 

[*£]„ = jf VB" ■ VBfrfu, , 

l<5r ( o] , = <*•■*> 

P?] ' = l V&, • VBpdw = {<*„, B"} v , 

[A= ! a h*' WJu =(+" B}/ ‘)« 

where the index ranges are k,£ = 1, • • ■ , M . . We note that the presence of the matrix tifg in 
the component of (2.14) and mass matrix is due to the form of the V- inner product. 

We point out that the Fourier-Galerkin technique leading to the matrix system (2.14) is quite 
general in nature and can be used to approximate wave solutions in cases where the analytic- 
solutions described at the beginning of this section are unavailable or impractical to use. The 
general techniques described here are also easily adapted to different boundary conditions and 
more complex systems involving acoustic components. Test cases illustrating the basic method 
are given in the following examples and the extension of these techniques to a more complex 
coupled structural acoustic system is presented in Section 4. 


k% l [ # M {t) i r o 

0 d M {t) + G?{t) 

(2.14) 
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2.1 Example 1 

As a first test of this approximation method, we considered the forced wave* equation 

(ptt = c 2 Aj) + //(/, r, (9) , (r,J)6fU> 0 , 

V0 ■ f = 0 , (r,0) e r ,* > 0 , 

0(0, r, 6) = 0<((U\0) = 0 , (r,0) G ft 

on the circular domain ft of radius a — 1. The true solution 0(i,r, 0) - F 2 [cos( 27 n*)- l] 2 [sin(30) + 
cos(0)] was used to generate the forcing function g , and the problem was then discretized to yield 
the matrix system (2.11). The absolute and relative errors obtained by integrating this system 
to time 7 = 1 and comparing with the known true solution on a uniform mesh are recorded 
in I able 1. I ho true solution at that time is depicted in Figure 1 while lino plots comparing 
the approximate and true solutions along the central line ( = {(.r, //) : - 1 < x < 1 , y — 0} are 
given in Figure 2. 

1 he solution for this problem was chosen in a manner that allowed us to fix the Fourier 
limit and examine the convergence as the (Jalerkin limit was increased. The rapid convergence 
exhibited by the results in Table 1 is consistent with that expected from the radial Legendre 
basis. A slight decrease in the convergence rate is noted when N = 15, and we believe that this 
is due to error in time discretization rather than spatial approximation (a standard fourth/fifth 
order Runge-Kutta routine was used to integrate the system). In this and other examples 
that, we have examined, it can be noted that as long as the solution being approximated is 
sufficiently smooth at the origin, the approximation method performs well and accurate results 
can be obtained with a relatively small number of basis functions. For example, the plots in 
Figure 2 demonstrate that while small oscillations near the origin are present in the approximate 
solutions obtained with N = 3 and N = 6, they are gone with larger discretization limits and 
if fact, tin 1 approximate solution obtained with A = 9 is graphically indistinguishable from 
the true solution, for less smooth functions (e.g., functions that are only continuous and 
differentiable at the origin), some loss of accuracy does occur, although, in applications this 
has not been a, factor since these types of discontinuities generally do not occur in the physical 
systems whose dynamics we are simulating. 


M N 

2M 

j|0irw.r 0npp|| 

Untrue J&tipp 1 1 

||0rru<*l| 

3 3 

54 

1.0067 - 0 

.1357 - 0 

3 6 

96 

o 

1 

lO 

.1959 - 1 

3 9 

138 

.1269-2 

.1710-3 

3 12 

180 

.8567 - 1 

.1155-4 

3 15 

222 

.1912-4 

.2578 - 5 


Table 1. 


Absolute 1 and relative errors when approximating the wave solution 
— / 2 [cos(27rr) — l] 2 [sin(30) + cos(0)]. 
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2.2 Example 2 

A second means of testing the accuracy of the mass and stiffness matrices in the system (2.14) 
and hence checking the accuracy and convergence properties of the approximation method is 
by discretizing the eigenvalue problem 

A<f> + 7 2 <J> = 0 , (r,0) G ft , 

V<J> . f = 0 , (r,0)e r 

(see (2.2)) which arises when separating variables in the homogeneous wave equation. As 
discussed previously, the eigenvalues j 7nn are determined by solving for the zeros of the nonlinear 
system ~ Tn d ^ a ^ = 0 where J m is the m th Bessel function of the first kind. Several values obtained 
with a — 1 are summarized in Table 2 (see page 343 of [9] for details). 

lo compare these results with those obtained via the Fourier-Galerkin expansion, it is noted 
that under approximation in weak form, the eigenvalue problem (2.2) yields the generalized 
matrix eigenvalue problem 

K$i) M = ~) 2 (2.16) 

with c 2 — a — 1 in the mass and stiffness matrices and respectively (see (2.15) for the 
definition of these matrices). The approximate eigenvalues obtained by solving (2.1(i) are sum- 
marized in Tables 3 and 4 for the limit choices M — 6, N — 6 and M = 6, N — 12, respectively. 
With the first choice of limits, it becomes difficult to distinguish the higher order eigenvalues 
and these are omitted from the table. We first note that with M = 6, N = 12, very accurate ap- 
proximations are obtained with the largest relative 1 error (4.4 x 10 -4 ) occurring when m — n — 6. 
Moreover, we see that while the definition of the modified Fourier component m (see (2.12)), 
which incorporates the analytic decay of the Laplacian near the origin, changes at in = 5, this 
has not reduced the accuracy of the method. Hence this example further illustrates the efficiency 
and accuracy of this approximation method for the wave equation on a circular domain. 
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m = 0 

m = 1 

m — 2 

in = 3 

m = 4 

in = 5 

in = fi 

0 


1.8412 

3.0542 

4.2012 

5.3176 

6.4156 

7.5013 

1 

4.8317 

5.3314 

6.7061 

8.0152 

9.2824 

10.5199 

11.7349 

2 

7.0156 

8.5363 

9.9695 

1 1 .3459 

12.6819 

13.9872 

15.2682 

3 

10.1735 

1 1.7060 

13.1704 

14.5858 

15.9641 

17.3128 

18.6374 

4 

13.3237 

14.8636 

16.3475 

17.7887 

19.1960 

. 20.5755 

21.9317 

5 

16.4706 

18.0155 

19.5129 

20.9725 

22.4010 

23.8036 

25.1839 

6 

19.6159 

21.1644 

22.6716 

24.1449 

25.5898 

27.0103 

28.4098 

Table 2. Values of 7 mn 

obtained from the Bessel condition 

4 AjT> ) _ n 

dr — U * 

n 

in — 0 

m = 1 

in — 2 

in = 3 

in = 4 

in ~ 5 

in = 6 

IT 


1.8412 

3.0542 

4.2012 

5.3176 

6.4156 

7.5013 

l 

3.8317 

5.3314 

6.7061 

8.0152 

9.2824 

10.5199 

11.7349 

2 

7.0159 

8.5378 

9.9714 

11.3472 

12.6823 

13.9874 

15.2709 

4 

10.2668 

1 1.7350 

13.1858 

14.6280 

16.0386 

17.4007 

18.7540 

4 

13.7669 

15.2709 

17.2673 

18.3866 




Table 3. Values of 7 mn obtainer 

1 with Aj 

t = 6, N 

= 6 Fourier-Galerkin basis functions. 
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n 

771 — 0 

m — 1 

m — 2 

m — 3 

m — 4 

m = 5 

771 = 6 

0 


1.8412 

3.0542 

4.2012 

5.3176 

6.4156 

7.5013 

1 

3.8317 

5.3314 

6.7061 

8.0152 

9.2824 

10.5199 

11.7349 

2 

7.0156 

8.5363 

9.9695 

11.3459 

12.6819 

13.9872 

15.2682 

3 

10.1735 

11.7060 

13.1704 

14.5858 

15.9641 

17.3128 

18.6374 

4 

13.3237 

14.8636 

16.3475 

17.7888 

19.1960 

20.5755 

21.9317 

5 

16.4717 

18.0159 

19.5131 

20.9730 

22.4018 

23.8043 

25.1860 

6 

19.6250 

21.1860 

22.6929 

24.1559 

25.5937 

27.0152 

28.4223 


Table 4. Values of 7 „ m obtained with M = 6, N - 12 Fourier-Galerkin basis functions. 


3 The Plate Equation 

A second area in which one must commonly approximate the dynamics of PDE-based models 
posed on circular domains occurs when numerically simulating structural dynamics. While the 
assumption of small amplitude vibrations often leads to linear models, these models are often 
sufficiently complex so as to warrant special numerical techniques to capture the physics of the 
system. As a motivating example, we consider a thin circular plate having piezoceramic patches 
bonded in pairs to its surface (see Figure 3). When a voltage is applied to the patches, stresses 
are generated which can be used to invoke in-plane forces and/or bending moments in the 
underlying structure [8]. In this manner, the patches can be used to control plate vibrations [2] 
or acoustic sound pressure levels when the plate is an active component in a structural acoustic 
system [3, 4, 5]. For the analysis which follows, we assume that the edge of the plate is clamped 
since this quite closely approximates the condition found in several applications of interest. 
We point out, however, that the techniques which follow are easily adapted to other essential 
boundary conditions if warranted by the physical model. 


Piezoceramic Patches 



Figure 3. A thin circular plate with piezoeceramic patches bonded in pairs to its surface. 
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For a thin circular plate of radius a having clamped edges, the strong form of the equations 
modeling the transverse motion is 


Pphw tt + — 


(P Mr 2 ()M 


.2 


Or 
dw 


+ - 


r dr 


1 dM( 
r dr 


2 (PM 




2 dM 


rO 


1 d 2 M ( 


+ rr~rPr + — — = f(P r ^) > 


r d/’cW r 2 cW 


<90 2 


(3-1 


w(t.,a,0) = ——(t,,a,0) = 0 
dr 

iu{ 0 , r , 0 ) - ?e 0 ( r, 0 ) , «j t ( 0 , 0 ) = w;, ( ?•, 6 ) 

where p p and // denote the density and thickness of the plate and w is the transverse displace- 
ment. The general moments are given hy 


where 
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MrO = l)( 1 — u) 


d 2 w 1 dw\ ( 1 (Pu) 1 d 2 w\ 

drdo ~7 2 He) + Cl ^ ~ v ’\i'drd0dt 


are the internal plate moments (including discontinuous changes in /), // and cp due to the 
bonding of the patches to the plate see [8]), and 


(Mr)* = (Me)rc = 

i= 1 

are the applied moments generated by ,s pairs of patches. With E denoting the Young’s modulus, 
the parameters u, l) = an< ^ C D represent the Poisson ratio, flexural rigidity, and damping 

coefficient for the plate/patch structure. Here y t (?’, 0) denotes the characteristic function which 
has a value of 1 in the region covered by the i th patch and is 0 elsewhere. Moreover, Ui(t) is the 
voltage into the i tk patch and K t is a parameter which depends on the geometry, piezoceramic 
material properties and piezoelectric strain constants (see [8] for details). We point out that the 
piezoceramic material parameters /C M i = 1,* • • ,.s as well as the plate parameters p p , D,cp and 
v should be considered as unknown and in applications must be estimated using data fitting 
techniques analogous to those discussed in [7]. The piecewise nature of the material parameters 
and input moments is one motivation for approximating the problem in a weak or variational 
form. 

To motivate the development of a Fourier-Galerkin scheme for this problem, we proceed in 
a manner analogous to that used in the last section and consider first a simplified version of 
the problem where separation of variables can be employed. To this end, we analyze first the 
equation 

Pph—^T + DV 4 w = 0 . 

di. 
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(3.3) 



which models the transverse motion of an undamped thin circular plate devoid of patches and 
having constant stiffness and density. Because the manner of solution is similar to t hat discussed 
in the last section, we highlight only those details which are necessary for motivating the form 
of the approximate solution for this problem. 

As before, we separate the temporal and spatial variables, in this case taking w(t,i\0) = 
T(t)\Y(r,Q), to obtain 

X/ 4 W = 7 4 14 / 

dw ( ;u ) 

W(a, 8 ) = — -(«,0) = 0 
or 

and T n + u> 2 T = 0 , t > 0. The separation constant is 7 4 = ppl ^ where u> again denotes the 
circular frequency. A second separation of variables then yields solutions of the form 


W{r) = + B m I m (^r) + C m V m ( 7 r) + D m K m {yr) (3.5) 

where J m and Y m are Bessel functions of the first and second kind and I m (z) = i~ m J 7 n (iz) and 
I\ m ~ lirn p ^ m ~[I^ p (z) — I P (z)\/ sin(p 7 r) are modified Bessel functions of the first and second 
kinds, respectively. We note that for fixed m and z — > 0, these functions have the limiting 
forms 

*,711 -,711 

Jm{z) ~ 7Z f Im(z) ~ — 7 

2 m m\ 2 m m! > 

2 m (m — 1)! „ , . 2 m (ni- H! 

Cn(~) ~ T~ ^ km{z) 


7TX 


2; 


(see [ 1 ]) which immediately shows that the Bessel function and modified Bessel function of the 
second kind grow unboundedly at the origin. Hence C m and D rn are taken to be zero in order to 
avoid infinite deflections and stresses at r = 0. Furthermore, the condition W{a,0) = 0 yields 
the condition 

U A) 


A m = -B„ 


J m ( A) 


(3.7) 


(3.8) 


where A = 7 a. Finally, the boundary condition. ^r(a, 6 ) = 0 requires that A solve 

= 0 

in order that (3.5) yields a nontrivial solution to (3.4). By solving (3.8), one obtains the eigen- 
values A mn and by combining (3.5) and (3.7) with the Fourier coefficients, the eigenfunctions 
or modes are found to be 


¥mM) = 


Jr, 


''77171 1 


Jm (^mn ) 


'A r S 
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(•!•!>) 


a ) Im(Kin) V (l 
With suitable smoothness assumptions on the the initial conditions, it then can be shown that 


the solution to (3.3) has the form 

oo oo 

w{t,r, 0 )= ^2 £>,„„(*) 


Jr, 


J ,n ( A ) 

U a mn ) 


/„ 


^77171 ^ 


imO 


(3.10) 


771 — — OO 11 — 1 

where the coefficients w mn (t) depend on the initial conditions. Again, we note that the motiva- 
tion given here is formal and we refer the reader to any standard PDE tex for details concerning 
this derivation. 
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Having motivated the general form of the analytic solution for a specific case of the plate 
equation, we can now present a Fourier-Galerkin method for approximating plate dynamics in 
the more general cases. As noted at the beginning of the section, the weak form of the system 
equation (3.1) will be used to accommodate potentially discontinuous material parameters. 
Moreover, the use of the weak form accommodates unbounded input operators (as is the case 
when piezoceramic patches are used as control elements) as well as reduces the smoothness 
requirements on the basis functions. 

The state for the problem is taken to be w and the Hilbert space // — L 2 (F o) with the 
energy inner product 

( w ^})h = / Pv hwT 1 d 7 

Jr 0 


is used as the state space. We also define the Hilbert space of test functions V = //q(T o) = 
{0 G // 2 ( I q ) : p(a,0) — ?/>'(«, 0) — 0} with the inner product 


( uk 7 1 ) v = { Mr 





dh ) 


do 2 


+ 2 f — M r 0 




’ drdO J 


dy 

' do 


where cp = 0 in the moments appearing in this definition and (F, (!) = fy. F(id~j with 
ihj = vdrdO. As detailed in [6], a weak or variational form of the equation describing the motion 
of a damped thin circular plate having .s piezoceramic patch pairs is 

<h 

( 3 . 11 ) 

= f 'Y F.,v, t (l)\ t (r, 0)X" 2 ijd^ + f fT/d-f 

l== | ^Fo 

for all // G V. Again, the internal moments are given in (3.2). 

To obtain an appropriate form for the approximate solution, we begin with the Fourier 
expansion in 0 

w A (t,r,0)= Y Fu m (t,r)e ,mB . 

771 = — M 

As was the case when describing the wave approximation, we will use the complex Fourier 
expansion while describing the method since it simplifies the notation. Due to the potential 
presence of patches or other actuators on the plate, however, the method is more easily imple- 
mented using the real (trigonometric) expansion (with patches present, the complex expansion 
leads to a complex system matrix which proves troublesome when solving the Riccati equation 
in the control problem). We have omitted details concerning the real expansions used in the 
implementation since it is straightforward to interchange between the two expansions. 

To determine an appropriate expression for w m (/,r), we first note that it must satisfy the 
conditions 

w m (t,r) = 0 , m ^ 0 ; 2^ = 0 

Or 

at /• — 0 in order to guarantee uniqueness and differentiability (see the discussion about the 
analogous requirement for the wave solution). In light of these requirements, an appropriate' 
expansion for w m (t,r) is 

j\ rm 

iFn(i-r) = y «w(0' Jn,| «;; i (e 

n= 1 


I l, l ,l,M W ld 1 + j + j + / p 

■/] o ^ I o «/ 1 q # ^ 1 ^ I o 


i m , dij 


d 2 ij (Pij i h) 

1 r9 do 2 + h drdo 2 oo 
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where £?” l (r) is the n ik modified cubic spline satisfying £?™(fl) — dB ^ a ^ = 0 with the condition 

dB dr^ ~ 0 being enforced when m = 0 (this latter condition guarantees differentiability at 
the origin and implies that N m = N + 1 when m ^ 0 and N 7n = N when rn = 0, where N 
denotes the number of modified cubic splines). The total number of plate basis functions is 
Af = (2 M + l)(N + 1) — 1. The inclusion of the weighting term with 


m ~ 


0 

1 


771 == 0 
777 ^ 0 


is motivated by the asymptotic behavior of the Bessel functions (see (3.6)) as r — » 0. It also 
serves to ensure the uniqueness of the solution at the origin. The Fourier coefficient in the weight 
is truncated to control the conditioning of the mass and stiffness matrices (see examples in [(>]). 
To simplify notation, the approximate solution is written as 


M N™ 

r,«)= £ 

in — — M n=l 


M) 


A: = l 


(3.12) 


where Bg(r,6) = 7 -l ,n l B™(r)e tm6 . A comparison with (3.10) indicates that, as was the case with 
the wave approximation, the Bessel components of the analytic solution are essentially replaced 
by approximating elements suitable for the problem under consideration. In the wave problem, 
the approximate solution was expressed in terms of weighted Legendre polynomials whereas here 
the radial basis is comprised of weighted cubic splines. In both cases, the weights were used to 
incorporate the analytic solution behavior at the origin into the approximate solution, thus guar- 
anteeing that it had the correct continuity and differentiability characteristics at that point. 

To obtain a corresponding matrix system, the Af dimensional approximating subspace is 
taken to be H* = span{B g } and the product space for the first-order system is Ti x x 7T V . 
The restriction of the infinite-dimensional system (3.11) to the space 7i Af x 7"t A then yields the 
matrix equation 


• Kg 0 ' 

0 . 

■ Kg 0 - 
0 M* 




' 0 

yv^(t) 


. -Kg 

W*( 0 ) ' 


' sg ' 

VV^fO) 


. gg . 


Kg 


W' v (t) 


0 

tig 


u( t . ) + 


0 

Fg(t) 


where W^(t) = [wgt), w 2 (t), • • • , w^(t)] T denotes the Af x 1 vector containing the approximate 
state coefficients. The component matrices and vectors are given by 

Kg = K [)i 4- K[)2 + Kp 3 + h [)4 + A os 
K C p — Ii CD \ T A CD 2 T A cd 3 T A C d4 T A C d5 

(3.13) 

Kl, = («>»■<)„ • Kl, = (”'•<)„ 
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where 
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[ K Ds]f,k = 2 f D(\-u) 
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~w dl 


with the index ranges are kj = 1, • • • ,A/\ The matrices K cp i — A CjD 5 are defined similarly with 
the inclusion of the parameter cp in the various integrals. Finally, we remind the reader that 
P p , D,v and cp are piecewise constant in these definitions due to the presence of the patches. 

For application purpose, it useful to note that the matrix system for the plate can be written 
as the Cauchy system 


r (0 - A Ar y M (t.) + B*u{t) + G*(t) 
y>) -• vo 


(3.14) 


where y^(t) = [VV^f^VV^ ( t )] = [w\ (/), • • • , tcyvCO? (0> * * * ? ^yv(0] denotes the 2 A/' X 1 vector 
containing the generalized Fourier coefficients for the approximate displacement and velocity. 
In th is form, finite-dimensional parameter estimation and control problems can be readily dis- 
cussed. 


3.1 Example 3 

As an initial test of the Fourier-Galerkin approximation scheme, we considered the undamped 
steady state problem 


C ?Mr 2 dMr 1 dMe 2 O 2 M r0 2 BM r0 1 (PMp _ 

dr 2 r dr r dr r drdO r 2 dO r 2 d$ 2 

with the moments given by (3.2). In these expressions, the Poisson ratio, flexural rigidity and 
damping coefficient were taken to be v = |, D = | and cp - 0, respectively. The plate radius 
was taken to be a = .6 and the true solution iv(r, 6) = (cos(2irr/a) — l)sin(d) was used to 
generate the forcing function 


/(»', 0) = — 3 T 67rr/a + 167 t ,j r 3 /a 3 ) sin(27rr/a) 

+ (—3 + 127 r 2 r 2 /o 2 + 167r 4 ? ,4 /a 4 ) cos(27n-/a) 


sin (0) . 
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The absolute and relative errors for this example are summarized in Table 5. By choosing a true 
solution with a known Fourier coefficient, the Fourier limit could be fixed at M = 1 and the con- 
vergence rate could be observed as the Galerkin limit N was increased. For sufficiently smooth 
forcing functions, we would expect the method to exhibit 0(h 4 ) convergence as a result of the 
cubic spline basis [18, 19]. To check this, asymptotic errors were calculated by dividing the pre- 
vious relative errors by 16. Since the number of radial basis functions is doubled each time, this 
provides a means of checking whether or not the expected convergence rate is being maintained. 
By comparing the results in the 5 th and 6 </l columns of the table, see that 0(h 4 ) convergence 
is being exhibited by the method thus providing an initial test in the efficacy of the method. 


M 

N 

size(A^) 

ll^true W<ipp|| 

1 1 Wt rue ~ Wapp \ \ 

Asym. Error 

IKrueH 

1 

5 

15 x 15 

.4362 - 2 

.2181 - 2 


1 

10 

30 x 30 

.1454 - 2 

.7269 - 3 

.1363-3 

1 

20 

60 x 60 

.6146 - 4 

.3073 - 4 

.4543 - 4 

1 

40 

120 x 120 

.6717-6 

.3358 - 6 

.1921 - 5 


Table 5. Absolute and relative errors for approximate solutions for Example 1. 


3.2 Example 4 

A second means of testing the accuracy and efficiency of the approximation method is by dis- 
cretizing the eigenvalue problem (3.4) and comparing the approximate eigenvalues and eigen- 
functions with analytic values that have been calculated for a simplified structure. This also 
provides a means of determining analytic values of the natural frequencies of a uniform plate 
which can then be compared with values obtained numerically and experimentally for structures 
involving plates to which patches have been bonded [2] and plates that have been incorporated 
in structural acoustic systems (as described in Section 4). 

For a uniform undamped plate to which no patches are bonded, approximation of the 
eigenvalue problem (3.4) via the Fourier-Clalerkin method described in this section yields the 
generalized matrix eigenvalue problem 

= 1 4 M Af '0* r (3.15) 

with the mass and stiffness matrices and A ff defined in (3.13). The plate dimensions 
a = ,2286 m (9 in), p p = 2700 kg/m 3 , h = .00127m (.05 in), and parameter choices E = 
7.1 x 10‘° N/m 2 and v = .33 were used which then yields the flexural rigidity D = 13.6007 N • m. 
We point out that these choices are consistent with the dimensions and parameters of an 
experimental plate currently being used in the Acoustics Division, NASA Langley Research 
Center so that the frequencies determined here could be compared with those obtained obtained 
experimentally (again, see the application in the following section). 

By noting the relationship / = where / is the frequency expressed in hertz and u> is 
the circular frequency, the natural frequencies of the fixed circular plate can be written as 

r = _L/Tl 2 /Z. 

J 2n \a) V Pph 
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For the given dimensions and parameter values, several frequencies deriving from the Bessel 
solutions A to the nonlinear equation (3.8), as reported on page 8 of [13], are given in Table (i. 
Approximate frequencies obtained by solving the generalized matrix eigenvalue problem (3.15) 
with the basis limits M = 6 and N = 24 are tabulated in Table 7. We point out that in these 
tables, the Fourier number rn can be interpreted as the number of nodal diameters while n is 
the number of nodal circles, not including the boundary. In comparing the Bessel and Galerkin 
results in the two tables, it can be seen that the Fourier-Galerkin provides accurate mass and 
stiffness matrices which translates to accurate approximations to the natural frequencies for 
the plate. With the basis limits used here, the largest relative error for the frequencies shown 
here is 8.1 X 10~ 4 when m — 3 , n = 6. Moreover, we see that while the definition of in in the 
weighting term changes at m — 1, the method loses no accuracy at that point. 

To qualitatively test the ability of the Fourier-Galerkin method to accurately approximate 
the decay of the solution in neighborhoods of the origin for increasing m, we compared the 
analytic eigenfunctions given by (3.9) with those obtained via the Fourier-Galerkin method 
with M — 5 and N — 16. Representative results using the two techniques are plotted in 
Figures 4 and 5 with additional examples demonstrating the 3-D behavior, corresponding 2-D 
slices and error results given in [6]. Qualitatively, the shape of the eigenfunctions in the two 
sets can be seen to be graphically identical and, by comparing the results for the (m = 5, n = 0) 
mode, it can be seen that the approximation scheme is accurately capturing the behavior of the 
solution near the origin. In combination with the eigenvalues results listed in the tables, this 
demonstrates that the approximation method accurately captures the physics of the problem 
throughout the circular domain with no loss of accuracy resulting at the origin, in spite of 
the coordinate singularity. Moreover, as demonstrated by the results in the last example, the 
expected convergence rates are also maintained through the use of the weight r l m L 


n 

m — 0 

m — 1 

m = 2 

m = 3 

m = 4 

m = 5 

m = 6 

0 

61.96 

128.95 

211.56 

309.58 

422.56 

550.38 

692.75 

1 

241.23 

368.90 

513.02 

673.33 

849.82 

1042.07 

1249.92 

2 

540.46 

728.34 

932.93 

1 154.26 

1392.14 

1646.34 

1916.70 

3 

959.46 

1207.39 

1472.15 

1 753.95 

2052.63 

2367.90 


4 

1498.20 

1806.12 

2131.29 

2473.02 




5 

2156.69 

2524.45 

2909.31 

3311.57 




6 

2934.91 

3362.52 

3807.60 

4269.79 





Table 6. Natural frequencies deriving from the Bessel expansions (in hertz). 
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m — 0 

m = 1 

m — 2 

m — 3 

m = 4 

m — 5 

in = 6 

0 

61.96 

128.95 

211.55 

309.52 

422.56 

550.38 

692.75 

1 

241.23 

368.96 

513.04 

673.40 

849.83 

1042.08 

1249.93 

2 

540.46 

728.35 

932.98 

1154.31 

1392.18 

1646.40 

1916.79 

3 

959.50 

1207.41 

1472.31 

1754.16 

2052.83 

2368.17 


4 

1498.37 

1806.28 

2131.37 

2473.57 




5 

2157.22 

2525.17 

2910.49 

3313.05 




6 

2936.35 

3364.46 

3810.17 

4273.26 





Table 7. Natural frequencies obtained via the Fourier-Galerkin scheme with M — 6, N — 24. 
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(m = 0, n = 1 ) (m = 1, n = 1) 



(m = 2, n = 2) 



(m = 5, n = 0) 



Figure 4. The (0, !),(!, 1),(2,2) and (5,0) modes obtained with Bessel expansions. 


(m = 0, n = 1) 


(m = 1, n = 1) 



(m = 2, n = 2) 


(m = 5, n = 0) 



Figure 5. The (0, 1),(1, 1),(2,2) and (5,0) modes obtained via the Fourier-Galerkin scheme 
with M — 5 and =16. 
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4 Application — Structural Acoustic Noise Control 

In t. tic* previous two sections, Fourier-fJalerkin methods for isolated acoustic and structural 
systems have been developed and illustrated through a set of examples. In this section, a 
structural acoustic application will be used to demonstrate the manner though which these 
techniques can be combined to yield an effective approximation scheme for simulating the 
dynamics of more complex systems. When modeling the experimental setup which motivates 
this application, cart 1 was taken to include the full coupling between the structural dynamics 
and the enclosed acoustic field, internal damping in the structure, and the effects of actuators 
such as piezoceramic patches which an 1 bonded to the structure. The full inclusion of these 
components leads to a model which accurately capture's the physics of the system but is difficult 
to approximate' using standard modal te'olmiques since mode's for the coupled system! (including 
coupling, damping and actuator effects) must themselves be approximated since' they are' ne>1. 
known analytically. Using a hourier-(Jalerkin technique, however, the approximate* system can 
be directly constructed with several of the* components actually consisting e)f matrice's from the* 
isolated acoustic and structural systems. In addition to illustrating the* me'thod for the* specific 
se'tup de'seribed here', this deunonstrates the flexibility of the method for approximating the 
dynamics for general structural acoustic systems having circular or cylindrical geometries. 

4.1 The Structural Acoustic System 

The structural acoustic system described here models an experimental se'tup currently being 
used for validation experiments in the Acoustics Division, NASA Langleyy Research Center, and 
the geometry and physical specifications were ediosen se^ as to be consistent with that apparatus. 
Specifically, t he* experimental apparatus is modeled by a cylindrical domain U having length f 
and radius a as pictured in Figure* (). At one end of the cylinder is a clamped flexible' plate of 
thickness h which is assumed to have Kelvin-Voigt damping. Bonded to the plate are sectorial 
piezoceramic patches which are' placed in pairs and excited out-of-phase* se) that a bending 
moment is produced whe'n voltage* is appled. The' patches and glue layer are* assumed to have* 
thicknesses / and 7)^, respectively. 



Piezoceramic Patches 



Figure 6. (a) The cylindrical acoustic cavity; (b) The circular plate with patches. 
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As discussed in [3, 4], an appropriate linearized variational (energy) form of the coupled 
system equations modeling this setup is 


f [ /»/V<A + 1 

n c J n J 

f gZdu 

Q 


(4.1a) 

- J r p p hw tt Y)d-y + Jr M r ^dj -f 

f -Meip-d'f 4 - / 

J r 0 t or J 1 

1 .. dh) 

— M r8T^d~f 

- 0 r l 0O 2 

(4.16) 

+2 / M re dj 2 f 

Jr Q r orov J r 0 

-2 M ^d 7 
r A 06 


(4.16) 

- Jr PS ~ u’tt) d~f 



(4.1 r) 

-- J^ 0)Vh)d-y + j 

f fydi 
r 0 


(4-1 d) 


for all appropriate test functions £ and r; (this will be clarified below). The differentials here 
are dto — r dr cW dz and r /7 = rdrdO with the overbars again denoting complex conjugation. 
We now consider the components of the equations and compare with those found previously. 
Throughout this section, the subscripts p and c will be used to denote plate and wave (cavity) 
components, respectively. 


Acoustic: With <J> G L 2 (tt) again denoting the acoustic velocity potential, the wave dynamics 
are contained in (4.1a) where the test functions £ are elements in iP(il) (see Section 2 
for details concerning the Hilbert spaces for this problem). When approximating the 
dynamics of the system, the 2-D basis described in Section 2 is tensored with an axial 
Legendre basis to yield the approximate solution 


Pc Me 

4, M (t,r,e,z) = Y. £ 


N?' m 


p=0 m — — M c n=o 

p+|m|+n^O 


where again, 


and 


m , [m\ = 0 , • • • , 5 

5 , \m\ — 6 , * - - , M c 


f P\(r) — 1/3 , p = ?/i = 0,w=l f N c + 1 . p+|m|/0 

P*' m (r) =\ , N*' m = < 

l F n (r) •> otherwise [ N c , p = m = 0 

Here P n (r) and P v {z) are the n th and p th Legendre polynomials which have been mapped to 
the intervals ( 0 ,a) and (0,/), respectively. We remind the reader that the term P\(r)— 1/3 
when p — m — 0 , n = 1 results from enforcing the condition f Q <p M (t, r, 0, z)(kv = 0 so as 
to guarantee that the functions are suitable as a basis for the quotient space while P m \ 
incorporates the analytic decay at the origin with m truncated to control conditioning. 
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Structural: Equation (4.1b) contains the internal plate dynamics while the external contribu- 
tions due to the perturbing noise source / and the excitation of the patches are contained 
in (4. Id) (compart' with (3.11) in the last section). As was the case when the isolated 
plate was considered, w £ L 2 (V o) denotes the transverse displacement and the test func- 
tions i] are taken in //q(I o)* When discretizing the system, the plate displacement is 
approximated by 0) as defined in (3.12) with M p and N p used here to denote the 

Fourier and spline limits, respectively. 

Coupling: The coupling between the plate dynamics and interior acoustic field is incorporated 
by including the backpressure />/<£*(/, r, 0) from the interior field as a force on the plate 
and assuming the continuity of velocity condition 

( / . , 0 , 0 ) = — w t ( / . r, 0 ) 
i )z 

at the surface of the plate (this latter condition is often designated the momentum con- 
dition in the literature). In the weak form of the system equations, these conditions are 
manifested as the first and second terms, respectively, in (4.1c). 


With the components thus described, we are now in a position to form the matrix system 
which results when the system dynamics are approximated. The Af = (2M V + 1 )( A r ;> + 1 ) — 1 and 
M — (2A/ r + 1 ) ( Ay + 1 ){P C + 1) — 1 dimensional approximating plate and cavity subspaces are 
taken to be H p = span {Bf }f_, and ll M = span{ Bf 4 } , respectively, where and Bf 4 a re 

the i th plate and cavity bases described above. Defining V — Af + A4, the approximating state 
spacer is H v — l I M x lip and the* product space' for the first order system is — H p X II P . 

The restriction of the infinite-dimensional system (4.1) to 7i p x 7 i v then yields the matrix 
system 

M r !l V (t) = A v y v (t) + B r u(t) + F v (t) 

M v y v ( 0 ) = ill . 

Here y v (t) — W V {t), i?^(/), \V v (t)] T , with the components t) v (t) = [6^ (/), • • • ,^vf(0] 7 

and W v (t) = denotes the approximate state vector coefficients while 

u(t) = [ui(/), • • • , u s (/)] T contains the s patch input variables. The system matrices and vectors 
are 
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We now separate the matrices into those containing acoustic, structural and coupling compo- 
nents to better illustrate the connections with those developed for the isolated acoustic and 
structural problems in the previous sections. 


Structural: The Af x Af component matrices M* , Kp and are the mass, stiffness and 
damping matrices which arise when solving the damped plate equation with fixed bound- 
ary conditions while the Af x 1 vectors and ffw are the corresponding control and 
forcing terms (see (3.13) for the various definitions). 

Acoustic: The M x M matrices M$ z and K$ z are the mass and stiffness matrices which 
arise when solving the uncoupled wave equation with Neumann boundary conditions on 
a cylindrical domain. As a results of the tensor properties of the 3-D basis, they ('an be 
succinctly defined as the tensor products 

M$ s = Af ^ 05 M% , 

K$ s = Mf 4 oo K% + I<? 05 M% 

where and are given in (2.15) and the ( P c + 1 ) x ( P c + 1 ) matrices AT and K z 
are defined by 

= [pMPA-y- . \ka„ = [i’;w’M) d: ■ 

The construction of M^ z and K^ z is completed by updating the row and column affected 
by the alterations used to guarantee that the functions are a basis for the quotient space. 
Finally, is a forcing term which incorporates any acoustic sources (see (2.15)). 

Coupling: The contributions from the coupling terms are contained in the matrices 



where Bf and Bf 4 are plate and acoustic basis functions, respectively, and the index 
ranges are k,( = 1 , • • • , Af and i, p = 1 , * • • , A/\ 

Initial Conditions: The vector y £ = [g^ 4 ,< 7 ^ , </^] 7 ( see (2-15) and (3.13)) contains the 
projections of the initial values into the approximating finite dimensional subspaces. 

We point out that several of the component matrices are identical to those defined when 
considering the isolated acoustic and structural systems while others can be efficiently con- 
structed from those components through basic operations such as tensor products. In fact, 
only the coupling matrices and A^ 2 containing quadrature values for the plate and acoustic 
bases must be constructed solely for this problem. This potential for decomposition into exist- 
ing structural, acoustic and coupling matrices also exists for more complex structural acoustic 
systems when a Fourier-Galerkin approximation scheme is used and is a further advantage of 
the method in coupled problems involving linear or weakly nonlinear components. 

4.2 Example 5 

As a first test demonstrating the convergence of the method for the fully coupled system, we 
considered the problem (4.1) with the dimensions a = . 6 , h = .00127,^ = 1.1 and parameter 
choices pj — 1 . 21 , c = 343, p p = 2700, v — .33, D — V2 (\-u 2 ) ~ 1 3.6007 and cp = .0001 1222 (the 
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effects of the patches are neglected in this example but will be incorporated in Example b). 
The true solutions 

w(r, 0) = / 3 (cos(27n7rt) - 1 ) sin 0 

(p{t,r,6,z) = -£f 2 (cos(2?rr/a) - 1 )z{( - z) 2 s\n0 

were used to generate forcing functions / and g on the plate and in the cavity, respectively (see 
the convergence tests in Sections 2 and 3). The resulting matrix system (4.2) was then integrated 
to time T = 1 and the absolute and relative errors in the plate displacement, potential and 

pressure at that time were calculated. The results for four sets of CJalerkin limits are reported 
in Tables 8 9 and 10. In these tables, the subscripts p and c. in the tables again refer to the 
plate and wave (cavity) indices, respectively, and the measurements were made on a 30 X 30 
grid on the plate and a 10 x 10 x 10 grid in the cavity. As demonstrated by these results the 
method for the fullv coupled system is exhibiting a convergence rate similar to that noted when 
the individual components were tested. Moreover, sufficient accuracy is obtained with tractable 
matrix sizes so that physical details are captured when the approximation method is used to 
simulate system dynamics and develop feedback control techniques. Finally, we point out that 
in spite of the fact that we are approximating nonaxisymmetric solutions, the computations tor 
this example could be performed on a Sparc 10-class workstation. 


Table 8. 


N e s;kc( A v ) 


2 2 
4 4 

6 6 


86 x 86 
182 X 182 
356 x 356 


8 608 x 608 


.4366 - 5 
.4294 - 5 
.2674 - 6 
.1588- 7 


.2207 - 2 
.2171 - 2 
.1352 - 3 
.8030 - 5 


Absolute and Relative Errors in Plate Displacement at 7 — .1. 


M v N p 


/(. .S i Zf( A V ) j | f Papp | 


2 86 X 86 

4 182 X 182 

6 356 x 356 

8 608 x 608 


.3132-2 
.1545 - 3 
.9554 - 5 ■ 
.2251 - 6 


I ^Iruf ^qppll 
Untrue |l 

.3569 - 0 
.1760- 1 
.1089 - 2 
.2565 - 4 


Table 9. Absolute and Relative Errors in Potential at T - .1. 


M p N p M c N c 


■S ize(A p ) 


2 2 86 X 86 

4 4 182 x182 

6 6 356 x 356 

8 8 608 x 608 


.7573 - 1 
.3739 - 2 
.2312 - 3 
.5453 - 5 


[Ptrue ~~Papy>H 

||ptT-udl 

.3566 - 0 
.1761 - 1 
.1089 - 2 
.2568 - 4 


Table 10. Absolute and Relative* Errors in Pressure at T — .1. 
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4.3 Example 6 

In this final example, the numerical techniques just described are used to simulate the system 
dynamics of the experimental structure described at the beginning of this section. A common 
technique for determining the dynamics and natural frequencies of coupled structural acoustic 
systems is to subject the structure to an impulse (through an impact hammer hit or a voltage 
spike to the patch) and measure the resulting time and frequency responses. By exciting the 
system in this manner, a wide spectrum of frequencies can be excited and various system 
properties determined. Because there is no acoustic source inside the cavity, <, = 0 in (4. "la). 

The hammer impact can be numerically simulated by a triangular force /"to a point on the 
plate while the voltage spike can be approximated by a short duration triangular voltage u; in 
both cases, the approximate system response is calculated via (4.2). 

In order to remain consistent with the experimental setup being modeled, the length and 
radius of the cavity were taken to be 1.0668 m (42") and « = ,2286 m (9"), respectively with a 
plate having thickness h = ,00127 m (.05") mounted at one end. A pair of circular piezoceramie 
patches having thickness T = .0001778m (.007") and radius rad = ,01905m (.75") were located 
at the center of the plate (see Figure 7). The physical parameters that were chosen for the 
structure and acoustic cavity are summarized in Table 11. The flexural rigidity D for the plate 
was obtained using the “handbook” value E = 7.1 x 10 ,o Af/m 2 for the Young's modulus of 
aluminum The remaining choices are comparable to values found when estimating parameters 
tor the isolated plate with a similar patch configuration [2], We re-emphasize that in general, 

1 i^iS ra,netCTS mUSt be determined through parameter estimation techniques in order for 
the PDE model to fit the actual physical system under consideration. As demonstrated in M 
the Fourier-Clalerkin scheme presented here performed well when incorporated in a fit-to-data 
routine for determining physical plate parameters thus yielding a model which could be used 
ioi implementing model-based feedback control techniques. 




a = .2286 


Figure 7. Th 


acoustic, cavity with a pair of centered circular patches, the impact point 
/I - (.13, f) and the observation points pi = (.125, and cl = (. 125, If, .05). 
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Structure 

Acoustic Cavity 

Parameter 

Plato 

Plate + Pzt 

Parameter 

Cavity 

p p • Thickness (kg/m 2 ) 

3.429 

3.489 

PI {kg/m 3 ) 

1.21 

D ( N * m ) 

13.601 

13.901 



C[) ( N * m * see) 

1.150-4 

2.250-4 



V 

.33 

.32 



K h (N/V) 


.0267 






c (m /sec) 

343 


Table 11. Physical Parameters for the structure and acoustic cavity. 


Voltage Spike to a Centered Patch Pair - Axisymmetric Response 

As a first example demonstrating the response of the system model to a simulated impact, 
we applied a triangular voltage spike u(t) as input to the patches and integrated the matrix 
system to obtain a time history of the response-. Because the patches are centered on the 
plate*, the response was uniform in 0 which implied that the Feuirier limits could he* taken to 
be M v = M c = 0. The remaining limits were taken to be N p = 12, Ay = 9 and P . = 9 
(the'se value's resolve the* range* of frequencies being e*xamined). A time history of the system 
response throughout the interval [0, .5] was calculated at the plate* and cavity points pi = 
(.125,4*) and cl = (.125, 4f , .05) (see Figure 7), and the resulting trajectories and frequencies 
are* plotted in Figure* 8. Phis temporal interval was chosen sine*e it was sufficiently long so as 
to demonstrate the system dynamics but short enough so that the higher frequency response's 
were not completely lost. The* off-center observation points were chosen to demonstrate the 
gene-rality of the* methoel anel provide a basis for comparison with the nonaxisymmetric results 
in the next e'xample. Finally, the system fre*quene'ies are* summarized in Table 12 with the 
notatiem p anel c being used to de'signate* those fre*e|ue*ncie*s wliie'h are observed at the* plate* 
point, /d and cavity point cl, respectively. 

To eletermine the e*ffe*cts of coupling, internal plate damping, and the presence of the patches 
eni the system, it is illustrative to compare the system results listed in Table 12 with those of 
the* isolated components. In Section 3, the natural frequencies for an isolated and undamped 
plate were obtained. The natural frequencies for the isolated wave equation in a cylindrical 
cavity can be determined by separating variables in the 3-D wave equation having Neumann 
boundary conditions. This leads to a 3-L) Helmholz equation which, after an analysis similar 
to that presented in Section 2, yields the natural frequencies 


J mnp 



(4.3) 


where p = 1,2. • • • , m = 0, 1,2, • • • , n = 0, 1.2, •• • and A mn = % nn a are zeros of (2.6) (see 
[3, 6, 9] for details). For the previously mentioned dimensions, the frequencies of axisymmetric 
modes (m = 0) are listed along with those for the isolated and undamped plate in Table 13 
(see Tables 2 and 6 for a compilation of the zeros y, un and natural frequencies for the isolated 
plate, respectively). 



By comparing the values for the individual plate and cavity in Table 13, with the system 
harmonics in Table 12, it can be seen that although the frequencies agree quite closely, there 
are slight differences due to the fact that the system (4.1) involves not only coupling between 
the plate and cavity but also includes damping in the plate. Specifically, the system frequen- 
cies associated with strong plate responses tend to be slightly less than those of the uncoupled 
and undamped plate while those system frequencies associated with strong acoustic responses 
are slightly higher than the natural acoustic frequencies of an isolated hardwalled cavity hav- 
ing these dimensions. Hence, through the Fourier-Galerkin method, physical characteristics 
about the system are determined and differences between the coupled system response which 
includes damping and the isolated, undamped components are illustrated. This latter obser- 
vation also indicates the difficulties which would be encountered if one were to attempt to 
simulated the coupled system dynamics using modes for the isolates! systems. 


Natural System Frequencies 



59 

P3' 

164 

c 

915 

p 

240 

P3* 

324 

c 

929 

p><* 

537 

Pi f 

483 

c 

971 



c 

045 





c 

807 




Table 12. System frequencies obtained with M p = 0, N p = 12, M c — 0 , N c = 9 and P c = 9 
basis functions; p - frequencies observed at the plate point pi = (. 125 , 4 tt/ 3 ), 
c - frequencies observed at the cavity point cl = (.125, 47r/3, .05). System freq- 
uencies can be compared with the ordered frequencies of the isolated and undamped 
plate and cavity given in Table 13. 


Plate 

( fmn ) 

Wave 

f fmnp ) 

(0,0) 

62 

(0,0,1) 

161 

(0,1,0) 

915 

(0,1) 

241 

(0,0,2) 

322 

(0,1,1) 

929 

(0,2) 

540 

(0,0,3) 

482 

(0,1,2) 

970 

(0,3) 

959 

(0,0,4) 

643 





(0,0,5) 

804 





(0,0,6) 

965 




Table 13. Axisymmetric natural frequencies for the isolated and undamped plate* and cavity 
(in hertz). 
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0 3r 


0.2 



Time (sec) 



Figure 8. The plate and pressure responses to a voltage spike. 


Noncentered Hammer Impact - Non-axisymmetric Response 

As indicated previously, a second means of exciting the system is through a hammer impact 
which can be modeled by a short triangular input force. To demonstrate the approxima- 
tion of general nonaxisymmetric system dynamics, this force impulse was applied at the point 
/I = (.130, ), and a time history of the system response for the time interval [0, .5] was cal- 
culated at the plate and cavity points pi — (.125, and cl = (.125, .05), respectively (see 

Figure 7). The resulting acoustic pressure and plate acceleration are plotted along with the cor- 
responding frequency responses in Figure 9 (the plate acceleration models data that would be 
experimentally obtained with an accelerometer, [2], and more clearly demonstrates the higher 
frequencies than docs the plate displacement). To resolve the frequencies below 1000 heitz, 
the basis limits were taken to be M p = 4,N P = 12, M c = 2, Ay = 9 and P e = 9 which yielded 
1230 coefficients in the system (4.2). As demonstrated by this choice of limits, the number of 
Fourier coefficients in the plate and wave expansions ('an be chosen to differ in order to reduce' 
system sizes although t are must be taken when doing so since the system response 1 contains 
contributions from both the plate and wave components. 

The system frequencies are summarized in Table 11 where again, the notation }) and r are 
used to designate those frequencies which are observed at the plate and cavity point pi and cl, 
respectively. For comparison, corresponding natural frequencies for an isolated and undamped 
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plate (see Table 6) and an isolated acoustic cavity (see (4.3)) are compiled in Table 15. By 
comparing the results in these table, it can be noted that six system modes corresponding to 
those of the isolated plate are excited with four of these responses being observed at both the 
plate and cavity points (we emphasize that because these are truly system modes, the remaining 
frequencies of 240 hertz and 422 hertz can also be measured in the cavity at various points; the 
response is simply weak at the point cl). Similarly, nineteen system modes corresponding to 
those of the isolated acoustic cavity are excited by the nonsymmetric impact with the strongest 
response in the asymmetric system mode having a frequency of 460 hertz (this corresponds to 
the (1,0,0) mode for the isolated cavity). As was noted in the previous example demonstrating 
the symmetric excitation, the system frequencies associated with strong plate responses tend 
to be slightly less than those of the uncoupled and undamped plate while those system fre- 
quencies associated with strong acoustic responses are slightly higher than the natural acoustic 
frequencies of an isolated hardwalled cavity having commensurate dimensions. We reiterate 
that these differences are due to the coupling between the plate and cavity as well as damping 
in the plate. This again provides motivation for considering a general Fourier-Calerkin method 
of the type described here when approximating the dynamics of a complex coupled system. 
The effect of coupling and damping are implicitly included in this approach while they must be 
explicitly (and occasionally artificially) enforced when employing modal approaches using the 
eigenfunctions of the isolated components. 


Natural Frequencies of the Coupled Structural Acoustic System 



Table 14. System frequencies obtained with M v = 4,N P = 12, M c — 2, N c = 9 and P c — 9 
basis functions; p frequencies observed at the plate point pi = (. 125, 47 t /3), 
c ~ frequencies observed at the cavity point cl = (.125, 47r/3, .05). System fre- 
quencies can be compared with the ordered frequencies of the isolated and un- 
damped plate and cavity given in Table 15. 


Plate (f mn ) 

Wave (fmnp) 

(0,0) 

62 

(2,0) 

212 

(0,0,1) 

161 

(0,1,1) 

929 

(1,0,5) 

916 

(0,1) 

241 

(2,1) 

513 

(0,0,2) 

322 

(0,1,2) 

970 

(2,0,0) 

729 

(0,2) 

540 

(2,1) 

933 

(0,0,3) 

482 

(1,0,0) 

439 

(2,0,1) 

750 

(0,3) 

959 

(3,0) 

310 

(0,0,4) 

643 

(1,0,1) 

468 

(2,0,2) 

797 

(1,0) 

129 

(3,1) 

673 

(0,0,5) 

804 

(1,0,2) 

545 

(2,0,3) 

874 

(1,1) 

369 

(4,0) 

423 

(0,0,6) 

965 

(1,0,3) 

653 

(2,0,4) 

972 

(1,2) 

728 

(4,1) 

850 

(0,1.0) 

915 

(1,0,4) 

779 




Table 15. Nonaxisymmetric natural frequencies for the isolated and undamped plate and 
cavity (in hertz). 
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Plate Acceleration at f 125. 4p»/3 v 


Amplitude Spectrum of Displacement at i, l2S.4pi,''3j 



Time (sect Hertz 

Figure 9. Tlir |)la.t.( ¥ and pressure 1 responses to an impart at (.ld.Tr/2). 

5 Concluding Remarks and Control Applications 

In this paper, a (la.Ierkin medhod for approximat ing the 1 dynamics of linear or weakly nonlinear 
systems in circular geHunetrms is pivseMite'd. While tin* method is presented for linear acoustic, 
structural and structural acoustic prohlemis. tin* techniques arc* sufficiently general so as to lu* 
applicable in a large* variety of applications involving tlu* approximation of linear or weakly 
nonlinear operators on circular domains. 

To take advantage 1 of periodicity. Fourier expansions arc* first used to approximate the tan- 
gential behavior of the solution. The 1 crux of t he method then involves the 1 (oust rind ion of 
radial basis functions by combining traditional approximat ing funct ions (finite 1 elemeuits, piece- 
wise splines, spectral polynomials) with terms of the* form /*l m l where /// is a truncated Fourier 
coellicient and the* coordinate 1 singularity is assumed to be* at the 1 origin. This latter term incor- 
porates the 1 decay propeud ies of the 1 Hesse 1 ! or analyt ic solution near the 1 origin and is include'd 
so as to guarantev that the 1 approximate solution is unique and differentiable at the coordi- 
nate 1 singularity. The t runcation of the 1 Fourie*r coedlicient provides a nievins of ('out rolling t he 
condit ioning of the resulting matrix system. 

As demonst rat ed by sewral numeuical example's, the 1 rexsuiting F()urie*r-( laleu'kin t.e'rh 11 ii|iic' 
is accurate 1 , eflicient and sufficiently flexible so as to be* easilv extended to complex cempled 
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systems. By varying the choice of radial basis functions, the method can be tailored to a 
variety of applications and boundary conditions while the inclusion of t he decay component r !m{ 
guarantees that accuracy is maintained throughout the entire circular domain for sufficiently 
smooth system inputs. Hence the Legendre polynomials and modified cubic splines, used in 
this presentation to approximate acoustic and structural dynamics, can easily be replaced by 
other spectral functions, splines or finite elements as warranted by the situation. 

A final advantage of the method arises when the Fourier-CJalerkin techniques are used to 
compute feedback gains in optimal control problems. As discussed in [5], where the method was 
incorporated in a control scheme for reducing noise in a cylindrical coupled structural acoustic 
system, the optimal control problem for reducing a periodically driven state z in the state space 
Ft can be posed under approximation as that of determining a suitable u in a control space ( 
which minimizes 

■J P (u) = t £ {(Q v z v (t) t z p (t)) H v + (K p u{t)Mt))u} dt (5.1) 

subject to z v satisfying the matrix system of dimension V which results when the equations 
modeling the system dynamics are discretized. Here r is the period and Q p and Rf* are matrices 
which can be used to weight various components of the approximate state 1 and control. As noted 
in [5], a suitable choice for Q p , when an energy inner product is associated with the state* spare 
Ft, is a diagonal multiple of the mass matrix M v which results when a (Jalerkin method is used 
to discretize the weak form of the system equations. Since tin* mass matrix constructed in this 
manner is the identity with respect to the energy inner produc t, Hie choice Q v — VM v , V a 
diagonal weighting matrix, in (5.1) minimizes a weighted measure of the state energy. Henc e 1 a 
Fourier-CJalerkin method of the type described here can also be 1 advantageous when calculating 
feedback gains for experimental and numerical implementation of optimal control techniques 
to problems posed on circular geometries. 
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transform methods for the problems are not investigated. While the specific method is developed for the 2-D 
wave equations on a circular domain and the equation of transverse motion for a thin circular plate, examples 
demonstrating the extension of the techniques to a fully coupled structural acoustic system are used to illustrate the 
flexibility of the method when approximating the dynamics of more complex systems. 
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